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Abstract 

We consider the static and dynamic models of Cournot duopoly with tax evasion. 
In the dynamic model we introduce the time delay and we analyze the local stability 
of the stationary state. There is a critical value of the delay when the Hopf bifurcation 
occurs. 
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1 Introduction 



During the last decades revenues from indirect tax have become increasingly important 
in many economies. Substantial attention has been devoted to evasion of indirect taxes. It 
is well known that indirect tax evasion, especially evasion of VAT, may erode a substantial 
part of tax revenues [2] , [4] , [5] . 

In [3] a model with tax evasion is presented. The authors consider n firms which enter 
the market with a homogenous good. These firms have to pay an ad valorem sales tax, but 
may evade a certain amount of their tax duty. The aims of the firms are to maximize their 
profits. The equilibrium point is determined and an economic analysis is made. 

Based on [1], [3], [7], [8], [10], in our paper we present three economic models with tax 
evasion: the static model of Cournot duopoly with tax evasion in Section 2, the dynamic 
model of Cournot duopoly with tax evasion in Section 3 and the dynamic model with tax 
evasion and time delay in Section 4. 

In Section 2, in the static model the purpose of the firms is to maximize their profits. 
We determine the firms' outputs and the declared revenues which maximize the profits, as 
well as the conditions for the model's parameters in which the maxim profits are obtained. 
Using Maple 11, the variables orbits are displayed. 

In Section 3, the dynamic model describes the variation in time of the firms' outputs 
and the declared revenues. We study the local stability for the stationary state and the 
conditions under which it is asymptotically stable. 

In Section 4, we formulate a new dynamic model, based on the model from Section 3, 
in which the time delay is introduced. That means, the two firms do not enter the market 
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at the same time. One of them is the leader firm and the other is the follower firm. The 
second one knows the leader's output in the previous moment t — r, r > 0. 

Using classical methods [6] , [9] we investigate the local stability of the stationary state by 
analyzing the corresponding transcendental characteristic equation of the linearized system. 
By choosing the delay as a bifurcation parameter we show that this model exhibits a limit 
cycle. 

Finally numerical simulations, some conclusions and future research possibilities are of- 
fered. 



2 The static model of Cournot duopoly with tax evasion 

The static model of Cournot duopoly is described by an economic game, where two firms 
enter the market with a homogenous consumption product. The elements which describe 
the model are: the quantities which enter the market from the two firms Xi > 0, i = 1, 2; the 
declared revenues Z{, i = 1,2; the inverse demand function p : R + — > R + ( p is a derivable 
function with p' (x) < 0, lira p (x) = 0, limp (x) = bi, (a\ G M, b\ el); the penalty function 

x^a\ x^O 

F : M + — ► R + (F is a derivable function with F' (x) > 0, F" (x) > 0, F (0) = 0); the cost 
functions Cj : R+ — > R+ ( C, are derivable functions with C[ (x^ > 0, C" > 0, i = 1,2 ). 
The government levies an ad valorem tax on each firm's sales at the rate ti G (0, 1) and 
q G [0, 1] is the probability with which the tax evasion is detected. 

The true tax base of firm i is Xip (x\ + X2) ■ Firm i declares Zi < Xip (x\ + X2) as tax base 
to the tax authority. Accordingly, evaded revenues of firm i are given by Xip (x\ + x 2 ) — Zj_. 
With probability 1 — q tax evasion remains undetected and the tax bill of firm i amounts to 
tiZi. The tax authority detects tax evasion of firm i with probability q. In case of detection, 
firm i has to pay taxes on the full amount of revenues, Xip (x± + x^) , and, in addition, a 
penalty F{xip(x\ + x 2 ) — Zi). The penalty is increasing and convex in evaded revenues 
Xip(xi + x 2 ) — Zi. Moreover, it is assumed that F (0) = 0, namely law-abiding firms go 
unpunished. 

The profit functions of the two firms are: Pi : — ► R+, i = 1,2, given by: 

Pi = Pi(xi,X2,z 1 ,Z2) = {l-q)[xip(x 1 +x 2 )-C i (x i )-t 1 z i \ + 

q [(1 - ti) Xip (xi + x 2 ) - d (x^ - F (xip (x 1 + x 2 ) - Zi)] . (1) 

The first bracketed term in (1) equals the profit of firm i if evasion activities remain unde- 
tected. The second term in (1) represents the profit of firm i in case tax evasion is detected. 

The firm's aim is to maximize (1) with respect to output Xi and declared revenues z%. 
This aim represents a mathematical optimization problem which is given by: 

max Pi, i — 1,2. (2) 

{Xi,Zi} 

From the hypothesis about the functions p, F,Ci,i — 1, 2 , we have: 

Proposition 1 The solution of problem (2) is given by the solution of the following system: 
dPi 

= [I - qt 1 - qF' (xipixx + x 2 ) - z i )][p(x 1 + x 2 ) + Xip' (x ± + x 2 )} - d(xi) = 

dP 

— - = -(l-q)t 1 +qF'(xip(x 1 +x 2 )-Zi) = 0, i = 1,2. (3) 
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In what follows, wc will consider the penalty function F (x) = — stix 2 , s > 1, the cost 

functions C, (xi) — CiXi, c, > 0, i = 1, 2 and the price function p (a;) = — . 

x 

From (3) we can deduce: 
Proposition 2 If ; — ^~T c i < C2 < < H—^—L Clj t/j e n the solution of system (3) is 



given by : 



qs + q - 1 " ~ 1-9 
c 2 (1 - h) 



ci(l-ti) 



\2 ' 



(ci + c 2 ) 

C2 i - g 



2 _ 7 I ^2" 
(C\ + C 2 ) 



ci + c 2 



ci + c 2 



1-g 



(4) 



For the parameters c\ = 0.3, c 2 = 0.6, g = 0.12, t\ = 0.16, the variations of the variables 
Zi, Z2, and the profits Pi, P 2 are given in the following figures: 



Fig.2.Plot(s,z2(s)} 





3 The dynamic model of Cournot duopoly with tax eva- 
sion 

The dynamic model describes the variation in time of output Xi (t) , i — 1,2 taking into 

dPi 

account the marginal profits — — , i = 1,2. Assume that each agent adjusts its declared 

OXi 

dPi 

revenue Z{ (t), i = 1,2 proportionally to the marginal profits — — , i = 1,2. Then, the 
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dynamic model is given by the following differential system of equations: 
dPi 

Xi(t) = ki-^ = ki{[l- qti - qF' (xipfa +x 2 ) - Zi)]- 

[p (xi + x 2 ) + Xip (xi + x 2 )] - Ci (xi)}, (5) 

Zi(t) = h i -— !L =hi[-(l-q)ti+qF'(xip(xi+X2)-Zi)], i = 1,2 
ozi 

with the initial conditions Xj (0) = x i0 , Zi (0) = z iai i — 1,2 and hi > 0, ki > 0, i = 1, 2. 
For F (x) — -stix 2 , s > 1 and d (x^ — axi, Ci > 0, i = 1, 2 system (5) becomes : 

Xl (t) = ki{[l-qt 1 -qst 1 (x 1 (t)p(xi(t)+x 2 (t))-z 1 (t))}- 

■ [p {xi (t) + x 2 (t)) + Xl (t)p' ( Xl (t) + x 2 (t))] - a} 
x 2 (t) = k 2 {[l-qt 1 -qst 1 (x 2 (t)p(x 1 (t)+X 2 (t))-Z 2 (t))]- 

■ [p {xi (t) + x 2 (t)) + x 2 (t)p' (xi (t) + x 2 (i))] - c 2 } (6) 

«i(t) = fci[-(l-g)ti + g«ti(a;i(t)p(a;i(t) + X2(t))-«i(t))] 

z 2 (t) - /i2[-(l-g)*i + g«*i(a;2(t)p(a;i(t) + a;2(t))-a;2(t))] 
Xi (0) = x i0 , Zi(0) = ZiQ, i = l,2. 

System (6) has the stationary state (x*, x 2 , z ii z 2 ) given by (4) . 

Let Mi (t) = xi (t) - x*, w 2 (*) = x 2 (t) - x\, U3 (t) = zi (i) - 2;*, u 4 (t) = z 2 (t) - z 2 . 

By expanding (6) in a Taylor series around the stationary state (x*, x 2 , z*, z 2 ) and ne- 
glecting the terms of higher order than the first order, we have the following linear approxi- 
mation of system (6) : 

Mi (t) = fci{ai ui (t) + a Q1 u 2 (t) + a aiu 3 (t)} 
u 2 (t) = fc 2 {6i wi (t) + b Q1 u 2 (t) + &001U4 (t)} 

ua(t) = h±{c w ui (t) + c iu 2 (t) + C001U3 (*)} (7) 
Ui(t) = h 2 {d w ui (t) + d iu 2 (t) + dooiiii (t)} 

where: 

c 2 

aio = -Qsh- ! — 2 + (1 - h) (2p' (xl + x* 2 ) + x\p" {x\ +x 2 )), 

(1 - ti) 

aoi = -qshxlp' {xl + x* 2 ) + (1 - h) (p' (xl + x* 2 ) + x* lP '' (xl + x* 2 )) , 

l — ii 

gstici gstic 2 gsiici 

aooi — -; —, 0001 — — , cooi — «ooi — —qsti, cio — — , (o) 

1 — 1\ 1 — £1 1 — 1\ 

ho = -qst lX * 2 p' (xl + x* 2 ) + (1 - h) (p' (xl + x* 2 ) + x* 2 p" (xl +x*)), 

1 — 1\ 

c 2 

601 = -g-sti - 2 , 2 + (1 - h) (2p' (xl + x*) + x*p" (x* + x^)) , 
(1 - ti) 

coi = qst\x\p' (xj +x 2 ) , dio = qst x x* 2 p (x* + x 2 ) , d i = qsh - ° 2 - . 

1 — El 
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The characteristic equation associated to (7) is given by: 

A 4 + m 43 A 3 + to 42 A 2 + TO41A + m 40 = (9) 

where: 

m 43 = -fcidio - k 2 b i - /11C001 - h 2 d m 

m 42 = kik 2 a w boi + (kia w + k 2 b i) (hiCooi + h 2 d Q oi) - hikia oiC W - 

-k 2 h 2 b Q01 d 01 + h x h 2 c o l^ooi - kik 2 a m b w 
m 41 = kik 2 a w b i (fticooi - h 2 d 001 ) - k 2 hih 2 c 001 d 001 b m + 

+/iifciCi a oi (fc 2 &oi + h 2 d 00 i) - kik 2 h 2 b 0Q1 a idi + 

+kik 2 h 2 a 10 b oidoi + k 2 hih 2 b oiCooid a i + 

+kik 2 a ib w (h 2 d oi + /11C001) - fcife^iaooibiocoi (10) 
^40 = fcife^i^2(aio^oiCooi^ooi — aooi^oi c iorfooi + aooi^ooi c io^oi + 
+aoi^ooiCooi^io — aooi^ooiCoirfio — 010^001 cooi^oi + 
+aooi^ioCoi^ooi — aoi&iocooi^ooi)- 

A necessary and sufficient condition as equation (9) has all roots with negative real part 
is given by Routh-Hurwitz criterion: 

Z?i = m 43 > 0, D 2 = m 43 m 42 — m 41 > 0, D 3 = m 41 Z? 2 — m^m^ > 0, D 4 = m 41 D 3 > 0. (11) 

From the Routh-Hurwitz criterion we have: 

Proposition 3 The stationary state of system (6) is asymptotically stable if and only if 
conditions (11) hold. 

4 The dynamic model with tax evasion and time delay 

In [7] and [8] we have studied the rent seeking games with time delay and distributed 
delay. In the present section we analyze the rent seeking games with tax evasion and delay. 
For t ~ we obtain the model from [3]. For r = and t\ — we obtain the model from [1]. 
We consider the model from section 3 where we introduce the time delay r. We suppose the 
first firm is the leader and the second firm is the follower. The follower knows the quantity 
of the leader firm, X \ (t — t) , which entered the market at the moment t — r, r > 0. 
The differential system which describes this model is given by: 

xi (t) = k 1 {[l-qti-qsti{x 1 (t)p{xi{t)+x 2 (t))-zi(t))}- 

[p (xi (t) + x 2 (t)) + xi (t)p' (an (t) + x 2 (t))] - ci} 
X 2 {t) = k 2 {[l-qt 1 -qst 1 (x 2 (t)p(x 1 (t-T) + x 2 (t))-z 2 (t))}- 

[p ( Xl (t - t) + x 2 (t)) + x 2 (i) p { Xl (t-r)+ x 2 (i))] - c 2 } (12) 
Zl(t) = h 1 [-(l-q)t 1 +qst 1 (x 1 (t)p(x 1 (t)+x 2 (t))- Zl {t))} 
z 2 (t) = h 2 [-{l-q)t 1 +qst 1 (x 2 {t)p{x 1 (t)+x 2 (t))-z 2 (t))} 
x x (9) = ip(6), 9 e [-t, 0] , x 2 (0) = x 2Q , z t (0) = z i0 , h > 0, hi > 0, i = 1, 2. 
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For p(x) = — the stationary state of system (12) is given by (4) . 

With respect to the transformation u\ (t) = x\ (t) — x\, u 2 (t) = x 2 (t) — x 2 , u 3 (t) = 
z\ (t) — z\ , u 4 (t) — z 2 (t) — z 2 , and by expanding (12) in a Taylor series around the stationary 
state {x\,x 2 , z\, z 2 ) and neglecting the terms of higher order than the first order, we obtain 
the following linear approximation of system (12) : 



Ui(t) = ki{aioUx (t) + a 01 u 2 (t) + a oal u 3 (i)} 

u 2 (t) = k 2 {b w u 1 (t - t) + b iu 2 (t) + &ooi«4 (t)} 

U3{t) = hi{c w ui (t) + C01U2 {t) + C001U3 (*)} (13) 

Ui(t) = h 2 {d w ui (t) + d iu 2 (t) + dooiUi {t)} 

where a w , a i, a oi, &10, b i, 6001, C10, c i, c oi, dio, doi, dooi, are given by (8) . 
The corresponding characteristic equation of (13) is : 

A 4 + n 43 A 3 + n i2 X 2 + n 41 X + n m + e~ Ar (n 22 \ 2 + n 21 X + n 20 ) = 0, (14) 

where 

n 43 = TO43, n 22 = -kik 2 a ib w 

n 42 = kik 2 ai b i + (fciaio + k 2 b i) (/11C001 + h 2 d 001 ) - /iifciaooiCro - 

—k 2 h 2 booidoi + hih 2 cooidooi 
n 41 = kxfoaxoboi (hxcoox - h 2 d oi) - fahxfocooidooiboi + 

+hikic w a ai (k 2 b i + h 2 d oi) ~ kik 2 h 2 b oia idw + 

+kik 2 h 2 a w b oid i + k 2 hih 2 b oic oidoi 
n 4 o = ^lfe^i^aiofroiCooi^ooi — aooi&oicio^ooi + aooi^ooi c io^oi + 

+aoi6ooiCooi^io — oooi&ooicoio'io — oio&ooicooio'oi) 
n 2 \ = kik 2 a ibi (h 2 d oi + /iiCooi) - fcife/iiaooifriocoi 
n 2 o = k 1 k 2 h 1 h 2 (a oibioc id 00 i - aoi^iocooi^ooi)- 

The roots of (14) depend on r. Considering r as parameter, we determine r so that 
A = ilu is a root of (14) . Substituting A = iu> into equation (14) we obtain: 
u 4 — in 43 uj 3 — n 42 uo 2 + in 44 uj + n 4 Q + (—n 22 u) 2 + in 2 iOJ + n 2 o) (cosojt — isinuor) = 0. 

From the above equation we have: 

lo 8 + r 6 L0 6 + r 4 L0 4 + t 2 lo 2 + r = (15) 

where 

re = n>l 3 - 2n 42 , r 4 = n\ 2 + 2n 40 - 2n 43 n 41 - n\ 2 , 
r-i = n\ x - 2n 42 n 40 + 2n 22 n 20 - n\ x , r = n\\ - n\ {y 

If uj is a positive root of (15) then there is a Hopf bifurcation and the value of t is 
given by: 

1 0104^0 + 0203 

t = — arctq , (16) 

uj — aia 3 + a 2 a 4 iv 

where a\ = u 4 - n 42 uo 2 + n 40 , a 2 = — n 43 uj^ + n 41 uj , a 3 = n 22 u 2 - n 2 o, 04 = n 2 iWo- 
We can conclude with the following theorem: 
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Theorem 4 (i) If ojq is a positive root of (15) and Re(^)\ = i iA)a . T=Ta ^ 0, where tq is given 
by (16), then a Hopf bifurcation occurs at the stationary state (x\, x 2 , z\, z 2 ) as t passes 
through r . 

(ii) If conditions (11) hold and no > 0, then the stationary state is asymptotically stable 
for any t > 0. 

5 Numerical simulation 

For the numerical simulation we use Maple 11 and the following data: q = 0.3, s = 40, 
h = 0.16, ci = 0.2, c 2 = 2, fci = 0.05, k 2 = 0.01, h x = 0.05, h 2 = 0.01. The stationary state 
is: x\ = 0.34710, x* = 0.0347, z{ = 0.85075, z* = 0.03257. 

For t = the Routh-Hurwicz conditions are satisfied. Then, the stationary state is 
stable. 

The positive solution of (15) is lu = 0.010083 and r = 164.5979. For r G (0,t ) the 
stationary state is asymptotically stable and for re [t , oo) the stationary state is unstable. 
For r = To there is a Hopf bifurcation. 

For q = 0.3, s = 40, h = 0.16, Ci = 0.2, c 2 = 1.5, k x = 0.05, fc 2 = 0.01, h x = 0.05, ft 2 = 
0.01. The stationary state is: x\ = 0.4359, x* 2 = 0.05813, z\ = 0.824019, z* 2 = 0.059313. 

For t = the Routh-Hurwicz conditions are satisfied. Then, the stationary state is 
stable. 

The equation (15) has no positive root. Then, the stationary state is asymptotically 
stable for any r > 0. 

6 Conclusions 

In the static model with tax evasion, the parameters q and s characterize the behavior 
of the firms with respect to evasion. The presented figures allow the analysis of the declared 
revenues and the profits with respect to s. 

For the dynamic model with tax evasion, using Routh-Hurwitz criterion we have deter- 
mined the conditions for which the stationary state is asymptotically stable. 

For the dynamic model with tax evasion and time delay, using the delay r as a bifurcation 
parameter we have shown that a Hopf bifurcation occurs when r passes through a critical 
value ro- 

The direction of the Hopf bifurcation, the stability and the period of the bifurcating 
periodic solutions will be analyzed in a future paper. 

The findings of the present paper can be extended in an oligopoly case. 
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